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We consider a class of states in an ensemble of two-level atoms: a superposition of two distinct 
atomic coherent states, which can be regarded as atomic analogues of the states usually called 
Schrodinger cat states in quantum optics. According to the relation of the constituents we define 
polar and nonpolar cat states. The properties of these are investigated by the aid of the spherical 
0^ \ Wigner function. We show that nonpolar cat states generally exhibit squeezing, the measure of which 

0^ ■ depends on the separation of the components of the cat, and also on the number of the constituent 

atoms. By solving the master equation for the polar cat state embedded in an external environment, 
we determine the characteristic times of decoherence, dissipation and also the characteristic time of 
5— ( ■ a new parameter, the non-classicality of the state. This latter one is introduced by the help of the 

Wigner function, which is used also to visualize the process. The dependence of the characteristic 
, times on the number of atoms of the cat and on the temperature of the environment shows that the 

decoherence of polar cat states is surprisingly slow. 

PACS number(s): 42.50.-p, 42. 50. Fx, 03.65.Bz 
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The question why macroscopic superpositions are not observable in everyday life has been raised most strikingly 
by Schrodinger in his famous cat paradox. Recent experiments however, show that at least mesoscopic superpo- 
I sitions can be observed in quantum-optical systems. In quantum optics one usually speaks of a Schrodinger cat (SC) 
' state if one has a superposition of two different coherent states of a harmonic oscillator. In one of the experiments 
101 a superposition of two different coherent states have been created for an ion oscillating in a harmonic potential. 
In the other one ||^ two coherent states of a cavity mode were superposed, and also the process of the decoherence 
, between these states could be followed by monitoring the field with resonant atoms. The unusual properties of such 
■ states have been discussed theoretically in several publications, see e.g. [^|-||. 
^ ' A different type of Schrodinger cat like state can be created in principle in a collection of two- level atoms, as first 
O^' proposed in The terminology we may use, is the following: The individual two- level atoms can be regarded as 
^ ' the "cells" of the cat, and the cat is definitely alive, if all of its cells are alive, i.e. they are in the |-|-) state, and it is 
•I— I , definitely dead, if all the cells are in the ill, |— ) state. In the case of N atoms a prototype of a SC like state is then: 

ivi'sc)-^(i+, +,...,+) + !-,-,. ..,-)), (1) 

where each of the terms contain N pluses and N minuses. We shall call this state the polar cat state, because the 
two components are in the farthest possible distance from each other. This state is in the totally symmetric iV -I- 1 
dimensional subspace of the whole 2^ dimensional Hilbert space, and if such states are manipulated by a resonant 
electromagnetic field mode with dipole interaction, then the atomic system will remain in this subspace. This is the 
arena of the collective interaction of the atoms and the electromagnetic field, called superradiance [sj-p^. In this work 
we present results concerning the properties and dynamics of polar cat states (|l[), and also of more general collective 
atomic states, the generation of which have also been considered recently | |ll] , [l^ . 

Our approach of discussing the properties of quantum states like l^'sc) is based mainly on the method of the Wigner 
function, which is one of the possible quasi-probability distributions. It has become a customary tool for investigating 
quantum states of an electromagnetic mode oscillator, or an ion oscillating in an appropriate trapping field p3| , |l^ . 
The method of Wigner function is much less exploited, however, in the description of atomic states like (|l|). That is 
why we first summarize the essentials of this method, and then turn to the determination of the Wigner function for 
the cat state (|l|) in Section II. Next, in Section III. we consider more general cat like states, which we call "nonpolar 
cats" and deteremine their squeezing properties. Finally, in Section IV. we write down and solve the master equation 
for a cat state in an environment with finite temperature. We define and determine the dissipation and decoherence 
times of the system, and the characteristic time when the system becomes essentially classical. 
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II. THE WIGNER FUNCTION OF THE POLAR CAT STATE 



The TV-atom dipole interaction with the electromagnetic field is equivalent to the dynamics of a spin of j — N/2, and 
the phase space of the atomic subsystem is the surface of a sphere of radius \/ j{j + 1), (fi = 1), which is sometimes 
called the Bloch sphere. This phase space and quasiprobability distributions corresponding to various operators acting 
in the 2j + 1 dimensional Hilbert space have been introduced first by Stratonovich [l^ . Similar constructions have 
been considered independently by several authors p6|-p0|. We use here the construction and notation introduced 



by Agarwal ||17[. Similarly to the case of oscillator quasidistributions [£1 22|, the quasiprobability functions for 
angular momentum states are not unique either. Beyond the natural requirements that the possible quasiprobability 
distribution functions have to satisfy, there is a special property, called the product rule, that distinguishes the most 
natural choice among the possible quasiprobability distributions. This rule requires that the expectation value of a 
product of two operators could be calculated by integrating the product of the corresponding quasiprobabilities. This 
choice is essentially unique, and in accordance with most authors we call it the Wigner function for spin j. We note 
that the construction can be extended to include several values of j , and in the same spirit Wigner functions 

can be defined for arbitrary Lie groups |^^. We also note that it is possible to define joint Wigner functions for 
atom-field interactions, and then a fully phase space description of atom-field dynamics can be considered p6[ . Here 
we restrict ourselves to the problem of angular momentum with a fixed value of j. 

Using the procedure proposed in [ pT[ we shortly summarize here the method of quasiprobability functions in the 
2j + 1 dimensional Hilbert space. One first chooses an operator basis in this space, and the most straightforward 
set of operators is the set of the spherical tensor operators Tkq which transform among others irreducibly under the 
action of the rotation operators Their explicit expression is: 

f j K j \ 

where ^ ^ is the Wigner 3j symbol. They form a basis in the sense that any operator of the 



-m m 

Hilbert space can be expanded in terms of them and they fulfil the Hilbert-Schmidt orthonormality condition 
Tr {tI^qTk'Q') = Skk-Sqq'. 

Introducing the characteristic matrix of the density operator p with respect of this operator basis as: 

gKQ^T:^{pT]^Q) , (3) 
the Wigner function of the state p is defined as: 

2j K 



^p(^'<^) = V^^E E QKQYKQ{e,<l>). (4) 
K=OQ=-K 

The factor in front of the sum ensures normalization. We note that in a similar way one can associate a Wigner 

function Wa{0, (f>) to any operator A, by introducing its characteristic matrix: Axq = Tr (^ATj^Q^ , and then forming 

the sum as in Eq. (|^). It can be easily seen, that this is a very similar procedure according to which one introduces the 
quasidistributions of oscillator states and operators by the help of characteristic functions of the translation operator 
basis: D{a) = exp{aa^ — a* a) ||2l|,^ . The construction of Eq. (||) can be shown to satisfy the product rule mentioned 
above, giving the following result for the expectation value of an operator A: 



Tr(pA) ^]I^J^ J Wp{e,^)WAie,(l))smed9dcl). (5) 
For other types of quasidistributions of angular momentum, like the analogs of the oscillator P and Q functions see 

Similarly to the case of the oscillator, the Wigner function allows one to visualize the properties of the state in 
question. In the work of Dowling & al. |2^ graphical representations of the Wigner function of the number, coherent 
and squeezed atomic states were presented. The Wigner function of a cat state like (|lj) has been considered first in 

The characteristic matrix of the state given by (|l|) can now be calculated according to the definition, Eq. (^), taking 
into account that the density operator corresponding to |vl/sc) is 
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in the standard basis, with j = N/2. The characteristic matrix has the form: 



{qsc)k.q 



V2K + 1 



j K j 
-J 3 



2 w -j Q ] ' ^' ' V -j 2j -j 

and from Eq. (^) one obtains the foUowing resuh for the Wigner function: 

y^^{N-l)\{N + 1 + 1)1 



[YioiO)+Yioi7T-e)] 



fe, + (-l)^^Q,-2,) 



(2iV+ 1)! (sin6l)^cos(iV(/.) | 



(6) 



(7) 



(8) 



The first term, containing the sums of two spherical harmonics, corresponds to the individual states |+,+, . . . 
and | — ), while the last term arises from the interference term between the "living" and "dead" parts of (|1|) 
(the last two terms of the density operator) . 

Fig. 1 shows the polar diagram of this Wigner function for = 5 atoms. The two bumps to the "north" and 
"south" correspond to the quasiclassical coherent constituents, while the ripples along the equator - where the function 
takes periodically positive and negative values - are the result of interference between the two kets of Eq. (|l|). The 
factor cos{N(f)) in (0) shows that the number of negative "wings" along the equator is equal to the number of atoms. 
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FIG. 1. Wigner function for a polar cat state, Eq. (|8|) for the case of A'^ = 5 atoms. The absolute value of the function is 
measured along the radius in the direction {6,4>), and the surface is shown in light where the function is positive and in dark 
where it takes on negative values. 



III. NONPOLAR CAT STATES, AND THEIR SQUEEZING PROPERTIES 



A. Nonpolar cat states 



One can also construct more general SC states by taking the superposition of any two atomic coherent states. 
An atomic coherent state (a quasiclassical state) |l6j, |t) is an eigenstate with the highest eigenvalue m = j oi the 
component of the angular momentum operator 'pointing in the direction n': 
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(J.n)|r)=i|r). (9) 

The notation r refers to a specific parametrization of the unit vector n by its stereographic projection to the complex 
plane. It is connected with the polar angle /3 and the azimuth a of the direction n as r = tan(/3/2)e~*". The atomic 
coherent state can be expanded in terms of the eigenstates |j, m) of Jz p^ : 



2j 



1/2 ^j+m 



= E f-!M ^ sin^"+"(/3/2)cos^-"(/3/2)e-'(-'"+")"|j,m). (10) 

The superposition of two quasiclassical coherent states is given by the ket: 

Ivp,,) = ^l±h^. (11) 
V2(l + Re(ri|r2)) 

Recently Agarwal, Puri and Singh and Gerry and Grobe have proposed methods to generate such states in 
a cavity, via a dispersive interaction with the cavity mode. 



Z 




We choose here n = tan(/3/2), T2 = — n. Then /3 is the polar angle of the classical Bloch vector corresponding to 
the atomic coherent state |ti) {f3 is measured from the south pole), see Fig. 2. This means that the x component 
of the expectation value of the dipole moment in these states is proportional to ±(iV/2) sin/3, respectively, and the 
y component is zero. Any other equal weight superposition of two atomic coherent states can be obtained from this 
special choice by an appropriate rotation. The polar cat state of the previous section corresponds to the special case 
when the two points are the northern and southern poles of the Bloch sphere. If the centres of the two coherent states 
in question are not in opposite points of the sphere, then we will call their superposition as "nonpolar" cat states. 

The corresponding quasiprobability distribution functions still can be explicitly calculated. For the Wigner function 
of the cat state 1^12) one gets the following expression: 

w(e 6) = /^±1 V V ^^^^(^^')' V (-ip'^'"' + (-1)^^+'" + (-1)^^- + (-i)^^-g 

V 47r A.^^Z.^2(l + (cos/3)2.) ^(j + m)!(j - m)!(j + Q + m)!(j - Q - m)\ 

^ ( -rn- Q Q m) (sin/3/2)^(-'+")+'?(cos/3/2)20-")-Q yKQ(0,0). (12) 
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We present polar plots of this Wigner function in Fig. 3. for N = 5 atoms and for several values of 





FIG. 3. Wigner functions for the state 1*1/12) for N — 5 atoms, and for several different values of /3: (a) /3 = 20°, (b) /3 — 45°, 
(c) f} — 55°, (d) /3 = 70°. For smaller values of /3 the state goes over into a single coherent state, and then it has essentially 
only one positive lobe. This graphical presentation shows qualitatively that the y component of the dipole moment is squeezed, 
the maximal value of the squeezing in the present case (A^ — 5) comes about (3 — 43°. 



For small (3 values, the interference is weak and the maximum of the Wigner function is around 9 = 0. For larger 
f3-s the function has two maxima around 6 = ±/3, and the interference gets more pronounced. When (3 = n/2, the 
two maxima corresponding to the individual coherent states point in the x and —x directions, respectively. In this 
case we get back the Wigner function of the SC state of Eq. (^), rotated around the y axis by n/2. 
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B. Squeezing properties 



The expectation values of the dipole operators Jx and Jy are zero in the state ( |ll| ) with ti — tan(/3/2), T2 = — ti, 
which is a consequence of the mirror symmetry of this state with respect of both the x — z and the y — z planes. As 
it is known, the variances of the dipole operators, Jx and Jy are equal to each other in an atomic coherent state: 

{^Jx)\) = {^Jy)\)^J/2. (13) 

In order to calculate the variances in the state ( pi] ) , one can use directly expansion ( pO| ) and the known matrix elements 
of Jx and Jy, but the summations that occur are rather cumbersome to evaluate. A more effective procedure is to 
apply the method of generating functions ||l6|. All the necessary expectation values in a cat state can be calculated 
by the formula : 



where 



= {T,\JlJlJl\r2) 



^2 



{(1 



{(i + |nP)(i + |r2P)p- 



2i 



(AJ.)^ = ^ 1 



(2j_ 

1 + (C0S/?)2J 



(14) 



(15) 



is the (antinormally ordered) generating function. 

Inserting the necessary operators, we obtain the following expressions for the variances in the state given by (|Tl|): 

1) sin^ (3 



(16) 



(AJy 



(2j - l)(cos/3)2J- 



(17) 



1 + (cos/3)2J 

Comparing these results with Eq. (|l3|), we see, that except for some special cases the Jy quadrature is squeezed while 
the Jx quadrature is stretched in this state. The reason of this asymmetry lies in the fact, of course, that in the 
superposition ( pi] ) we have chosen states that are both centered in points lying in the x ~ z plane. 

One of the exceptional cases that is not squeezed is, if there is only one atom: j — 1/2. As it is easily seen, for 
j = 1/2 any state in the two-dimensional Hilbert space is a coherent state, and therefore it does not show squeezing. 
The two other exceptions are /3 = for any j, because then the two coherent states coincide, and (3 = 7r/2, which is 
the rotated version of the polar cat state. 

Writing Eq.(17) in the form (AJ^)^ = j(l — 5)/2, we can define the quantity S as the measure of squeezing. Analysis 
shows that if N is large enough, then the maximum value of S is 0.56 and it is achieved around (3m = 1.6/-\/iV. Figure 
4 shows the dependence of 5 on /3 for several values of TV = 2j. 



. 6 -1000 100 




FIG. 4. The /3 dependence of the quantity S in (AJ^)^ = j(l — (S(/3, j))/2, for several values of A'^ = 2j. S can be considered 
as the measure of squeezing for a cat state consisting of two atomic coherent states separated by the central angle 2(3 on the 
Bloch sphere. 
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IV. DECOHERENCE AND DISSIPATION 



As we mentioned in the Introduction, there have already been reaUzable methods proposed for the experimental 
generation of atomic SC states in a collection of two- level atoms ||ll|,|l2|. However, such an atomic ensemble can never 
be perfectly isolated from the surrounding environment. Further, any observation of these states necessarily leads to 
the interaction of the atomic system with a measuring apparatus. In both of these cases the atomic system interacts 
with a system containing a large number of degrees of freedom. A possible and successful approach to this problem 
considers that the static environment continously influences the dynamics of the atomic subsystem, which besides 
exchanging energy with the environment loses the coherence of its quantum superpositions and evolves into a classical 
statistical mixture. 

In this section we investigate the decoherence and dissipation of the atomic Schrodinger cat states embedded in an 
environment with many degrees of freedom, by writing down the master equation for the reduced density operator of 
the atomic subsystem. We provide the solution for the polar cat states (0). 



A. Model and solution 



We couple our ensemble of two-level atoms to the environment which is supposed to be a multimode electromag- 
netic radiation with photon annihilation and creation operators Uk and a^. Then the interacting system can be 
described by the following well known model Hamiltonian which considers dipole interaction and uses the rotating 
wave approximation: 

H = iVaJz + h'^ujkalak + '^9k (a[.J_ + akJ+^ , (18) 

k k 

where is the transition frequency between the two atomic energy levels, the ujk denote the frequencies of the modes 
of the environment and gk are the coupling constants. If we suppose the environment to be in thermal equilibrium at 
temperature T, then the time evolution of the atomic subsystem is determined by a master equation for its reduced 
density operator p{t) [ |30| , |3l| |: 

= -| (M + 1) (J+J-Pit) + P{t)J+J- - '2J-p{t)J+) - I {n) {J-J+p{t) + p{t)J-J+ - 2J+p{t)J^) (19) 

which involves the usual Born-Markov approximation and is written in the interaction picture. Here {n) — 
(exp(?iWa/(fcBT)) - 1)~ is the mean number of photons in the environment and 7 ~ {g{LjJa)a{u)a))'^ denotes the 
damping rate, where a is the mode density of the environment. 

Eq. (|l^) can be obtained also in a somewhat different context, as described in [Q. Then one assumes the atomic 
subsystem to be placed in a resonant cavity with low quality mirrors causing the damping of the cavity mode at a 
rate k. Under certain reasonable assumptions one can get Eq. ( p^ with 7 = 2 g(uj^Y / k. 

From Eq. (^9|) one can easily deduce the following equations for the matrix elements of the density operator 
PmAt) = {3,m\p{t)\j,l): 

= -1 [ (n) (2j-(j + 1) - ™(m + 1) - l{l + 1)) + ((n) + 1) (2j(j + 1) - m{m - 1) - l{l - 1)) ] p^.i{t) 
+ 7(") V(.?(J + 1) - ™("'' - 1)) (j(i + 1) - 1(1 - 1)) Pm-u-i{t) 



+ 7 ((n) + 1) V(j(j + 1) - Mm + 1)) (jXj + 1) - + 1)) P.^+U+l{t)■ (20) 

Thus the time evolution of a particular density matrix clement is coupled only to the two neighbouring elements in 
the corresponding diagonal for {n) > 0, and only to the neighbour with larger index at zero temperature. 

In the case of a polar cat state (consisting of = 2j atoms), the elements of the density matrix have zero initial 
values except for p-j^-j, pjj, P-j,j and Pj^-j {— P*-~j.j)- This implies that the density matrix elements, except for 
those in the main diagonal and for p-j^j and Pj,-j, remain identically zero for any time. Setting 7=1 (i.e. the time 
unit is 1 /"/) the equations for the elements in the main diagonal of the density matrix are the following: 



dpm,m{t) 
dt 



= - [{n) UU + 1) - m(m + 1)) + ((n) + 1) + 1) - m(m ~ 1))] p™,™(t) 

+ (n) + 1) - m(™ - 1)) Prn-l,rn-lit) 

+ ((n) + 1) (j(j + l)-m(TO + l)) (t) (21) 
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with the initial values Pm.m{t = 0) = ■^{Sm.j + Sm.-j) (cf- equation (^). The dynamics of p-j.j is governed by the 
particularly simple equation: 

^^^ = -j(2(n> + l)p-«W, (22) 
yielding immediately the following solution with the initial value p_jj (0) — 1/2 corresponding to the polar cat state: 

P-«W-^exp(-j(2M + l)i). (23) 



As expected, the stationary solution of Eq. (21) is the Boltzmann distribution of the stationary values pr, 

nio,\ i-c^^{-nuj{k^T)) _ ((„)/((„) + 



Pm,m = cxp -(m + j)— — /o- . in rj^v, = 7 (24) 



l-exp(-(2j + l);itJa/(fcBT)) ((n) + i) (l - ((n)/((n) + 1))'^'+') 



Approximate analytical time dependent solutions of Eq. (|T]) can be found especially for the case of superradiance, 
when Pjj{0) = 1 in |3^, see also and references therein. For the initial conditions corresponding to the polar cat 
state, the time dependent solution of equations (^ at zero temperature ((n) = 0) can be obtained by the following 
recursive integration: 



ftj (0 = ^exp (-2jt), 



(25) 



t 

1 



Pm,m{t) = exp (-brnt) \ - Sjn-j + &m+l / exp (&„<') Pm+l,m+l(i') dt' | , -j <rn< j 



where bm = j{j + 1) — m(m — 1). These equations show rather explicitly, how does the initial excitation cascade down 
to the zero temperature stationary state. 

For non-zero temperatures ((n) > 0) we have solved equations (21) numerically. We are going to analyze the 
solutions in the next subsection. 



B. Characteristic times 



Figure ^ shows the time evolution of the relevant density matrix elements pm,m{t), m = —j, —j + 1, . . . ,j (solid 
lines) and p-jj{t) (dashed hue), in the case of initial polar cat states consisting of 5 and 50 atoms, for (n) = 0, 1, 10. 
The actual value of p-jj{t) characterizes the coherence of the corresponding state, since P-j,j and Pj-j (= P-jj) 



are the only nonzero matrix elements outside the main diagonal. Their exponential decay (cf. Eq. (23)) is the 
decoherence, shown by the dashed lines in the plots of Fig. |^. Thus it is reasonable to define the characteristic time 
of the decoherence by 

_ 2 

^"""^ ~ iV(2(n) + l)' ^^^^ 

implying p^jj{tdcc) ^ p-j,j{0)/e. 

In contrast to the simple time dependence of p-jj{t), the dynamics of the main diagonal elements Pm.m{i) depend 
on the actual value of (n) and N rather sensitively. The zero temperature cases. Fig. g (a) and (d), clearly show the 
initial excitation, contained in Pj.j{0) = 1/2, cascading down to /9_j__j(oo) = 1 as given by Eqs. (p5|). At nonzero 
temperatures ((n) > 0) the time evolution of the Prn.m{t)-s is more complicated because of the coupling to both 
neighbours, cf. Eq. (pO|). 

More information can be extracted from the time evolution of the Pm,m{t)-s by calculating the energy of the atomic 
subsystem as the function of time: 

j 

E{t)={uj^J^){t)^uJi,Ti{p{t)Jz)^huJa, ^ mpm,mit). (27) 
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(a): N=S, <n> = 0, r = 1.27 (b) : N=S, <:n> = l, r = 2. 48 (c) : N=S, <:n> = 10, r = 3.06 




FIG. 5. Plots of the density matrix elements pm,m{t), m — —j, + (solid lines) and p-jj{t) (dashed line) versus 

time (the time unit is I/7). Plots are given for N = 5 atoms: (a, b, c), and = 50 atoms: (d, e, f), for (n) — (zero 
temperature): (a, d), for (n) = 1: (b, e), and for (n) = 10: (c, f). The solid lines may be identified as follows: in (a) the 
smaller m is the later the corresponding pm,m has its maximal value; in (b) and (c) the stationary values of the pm.m-s follow 
the Boltzmann distribution, see Eq. (p4|). In (d, e, f) the Pm,m-s follow each other along the axis labeled by the m. The dotted 
lines (starting from 1 at t = 0) show the normalized energy of the atomic subsystem 1 + E{t)/{—jhuJa.), see Eq. (p7|). The 
number r — tdiss/idcc is the ratio of the characteristic time of dissipation to the characteristic time of decoherence. 



The process of dissipation (i.e. the change of the energy of the atomic subsystem in time) can be very easily followed 
by studying E{t). This function, normahzed to the zero temperature stationary energy and shifted to vary from 1 
to its stationary value: 1 + E{t) / [—jfiujs), is shown in the plots of Fig. || by the dotted lines. Since its asymptotic 
behavior is exponential-like, it is reasonable to define the characteristic time of dissipation idiss by requiring 

\E{U,ss) - E{^)\ = \E{Q) - E{^)\/e. (28) 

In order to ensure that E{t) achieves its stationary value with a good accuracy in the plots of Fig. ^ we have set the 
time range to 5idiss- It is seen that the value of r = tdiss/^doc grows with both the temperature and the number of 
atoms. A more detailed analysis of this question follows later in this section. 

The initial state of the process, the polar cat state, has sharply non-classical features. On the other hand, at 
non-zero temperature the final stationary state of the present model is a thermal state, which is classical in its nature. 
(At zero temperature the stationary state is also non-classical, since it is the state |j, — j).) It is natural to ask, when 
does the transition from the non-classical to the classical stage occur? What is a good measure of non-classicality 
reflecting the change of non-classical nature of the corresponding state? 

The spherical Wigner function (^) provides a good answer to both of these questions. Quantum states are generally 
considered essentially non-classical if the corresponding Wigner function takes on also negative values. Therefore 
to answer the second question, for the measure of the degree of non-classicality we propose to use the quantity 
V ~ 1— (/+ — I-)/ ( J+ -I- /- ) , where /+ is the integral of the Wigner function over those domains where it is positive 
while /_ is the absolute value of the integral of the Wigner function over the domains where it is negative. Since the 
integral of the Wigner function over the sphere is 1, = 1, thus h' = 2/_/(2/_ + 1), and it is easy to see that 
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< v < 1. According to this definition, the bigger is the value of v, the more non-classical is the state, and for all 
classical states one has i/ = 0. 

Regarding now the first question, namely for how long is the state of the atomic system non-classical, we introduce 
a third kind of characteristic time inci- We define t^ci to be the time instant when the corresponding spherical 
Wigner function becomes non-negative everywhere on the sphere, i.e. becomes 0. Wc will return to this question 
in connection with the time evolution of the Wigner function, which we will present in the next subsection in more 
detail. 

Based on the information provided by the three kinds of characteristic times, we consider here the dependence of 
the process on the number of atoms and on the temperature. In Figure ^ we plot idiss (dashed line), tdec (solid line) 
and tncl (dotted line) as the function of the number of constituent atoms of the polar cat N, for several temperatures, 
on a log-log scale. 




1 10 100 1000 

Number of atoms 



FIG. 6. Plots of the characteristic times of decoherence, t^cc (solid line), dissipation, tdiss (dashed line) and non-classicality, 
tncl (dotted line) versus the number of atoms, on a log-log scale. The uppermost solid and dashed line are for (n) — (there is 
no plot for tncl at zero temperature, since the state stays non-classical), while the subsequent groups of the three kinds of lines, 
one under the other, are for (n) — 1, 10, 100, respectively. 

It is seen that the characteristic time of decoherence tdn c is inversely proportional to the number of atoms (the 
straight solid lines in Fig. ^), according to the definition (p6|). Compared to this, the characteristic time of non- 
classicality tncl decreases less rapidly with increasing number of atoms. The characteristic time of dissipation tdiss 
first slightly increases at non-zero temperature then it achieves a maximum which depends on (n) and finally it 
decreases nearly inversely proportional to the number of atoms. The values of tdiss at different temperatures seem to 
converge slowly beyond a certain number of atoms. 

It seems however rather surprising that the ratio tdiss /tdcc is not as large as such a quantity is usually expected to 
be 11,11 : in the case of iV = 1000 it is 4.04 for (n) = 0, and it is still just around 350 for (n) = 100. (Note that 
(n) — 100 corresponds to a temperature of 250 K in the case of typical experiments [^.) The ratio tdiss/tdec seems 
not even to vary considerably with increasing N beyond the maximum of tdiss mentioned above. Thus the process of 
decoherence is extremely slow in the case of a polar cat state which is coupled to the environment by an interaction 
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leading to the master equation (|19|). 

Similar effects have already been reported for other physical systems earlier ||3^ . In a recent work Braun, Braun 
and Haake investigated the decoherence of an atomic SC state |ti) + |r2) based on Eq. (|l^) for zero temperature. 
By evaluating a certain quantity characterising the decoherence rate at the initial time, and applying a semiclassical 
procedure for finite times they concluded that for atomic SC states with T1T2 = 1 the decoherence slows down. 

Our initial state, the polar cat state fulfils the former condition. The results presented in Fig. ^ derive from the 
solution of the master equation for the whole process. They are in agreement with the statements of Ref. where 
the initial stage of the decoherence is analyzed for the case of zero temperture. 



C. Wigner functions 

We illustrate now the process of decoherence and dissipation using the spherical Wigner function (^. 




-0.2 



-0.2 




c. 2 




K 0,1 





-0 , 2 



-0 , 4 



FIG. 7. Polarplots of the temporal change of the Wigner function representing the decoherence and dissipation of the initial 
polar cat state (^), made of 5 atoms and shown in Fig. ^ at zero temperature. The dynamics of the corresponding density 
matrix elements is shown in Fig. ^(a). The time instants are the following (in units of I/7): (a) 0.1, (b) 0.2, (c) 0.4 {= tdcc), 
(d) 0.506 (=fdiBs), (e) 0.8, (f) 2.5. 

In order to obtain its time dependence we have to calculate first the characteristic matrix QK,Q{t) — Tr (^p{t)Tl^ q 
from the matrix elements Pm,iit) according to 
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gK,Qit) = V2KTT j2 (-1)'"™ ( _m Q m - Q ) (29) 

m——j ^ ^ 

From Eq. ( p9| ) it can be seen that only gi^.o = 0,1,..., A^) and qn.n = {—^)^{0n.-n)* are nonzero. This 
fact (which is due to the initial conditions specified by the polar cat state) ensures that the azimuthal dependence 
of the spherical Wigner function is determined only by the real part of the spherical harmonic Ym,n{0,4') which 
is proportional to cos(-/V0). Therefore the Wigner function keeps its initial azimuthal symmetry during the whole 
process. Further, since QN.N{t) = (— l)^pj,_j(i), the azimuthal modulation of the spherical Wigner function explicitly 
shows the degree of the coherence of the actual state. 

Figures ^ and || show the polar plots of the spherical Wigner function transforming its shape in time at {n) — 
and in) ~ 10, respectively. The initial state is a polar cat made of A'^ = 5 atoms and its Wigner function is shown in 
Fig. 0. 




FIG. 8. Polar plots of the temporal change of the Wigner function representing the decoherence and dissipation of the 
initial polar cat state (^, made of 5 atoms and shown in Fig. |^, for (n) = 10. The dynamics of the corresponding density 
matrix elements is shown in Fig. |^(c). The time instants are the following (in units of I/7): (a) 0.01, (b) 0.019 (= idee), (c) 
0.031 (=tnci), (d) 0.045, (e) 0.058 (=tdiBB), (f) 0.25. 

In Figures and || the following main characteristics of the process can be identified. The decoherence is shown by 
the decreasing and finally disappearing ripples along the equator. The vanishing of non-classicality, i.e. the decrease 
of the parameter can be easily recognized as the decrease of the negative (dark) wings. At nonzero temperature 
they disappear exactly at tnci, as shown in Fig. ^(c). The dissipation is represented by the approach of the initial 
upper and lower bumps to each other. At zero temperature (Fig. ^) the upper bump disappears and goes over to the 
lower one. This stationary shape of the Wigner function corresponds to the lowest coherent state ]],—]) For 
(ri) = 10, when the stationary energy is close to the initial energy, not only the upper bump moves downwards but 
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also the lower one lifts upwards. The stationary Wigner function has nearly a spherical symmetry, although its center 
is not in the origin. 

In agreement with Fig. |[ the plots of Figures ^ and || show, that the timescales of the decoherence and of the 
dissipation are very close to each other in the case of 5 atoms, for zero temperature they are practically the same. 
Then the spherical Wigner function exhibits considerable azimuthal modulation (ripples) also at idiss- 

We may come back finally to the question of finding the characteristic time of non-classicality t^ci- According to 
the arguments given after Eq. (^9|) it is sufficient to study the Wigner function within a 0-range of length 2tt/N, e.g. 
< (j> < 2tt/N, because it is invariant with respect of rotations by (j) = k^, (fc = 1, 2, . . . , N) i.e. it has Cjv symmetry 
at all times. Therefore the spherical Wigner function of a polar cat state, while subject to dissipation and decoherence, 
has its minimum value always at cj) — tt/N. Thus in order to calculate t^ch it is sufficient to follow the time evolution 
of the section W{9, (p = tt/N). Further, in connection with the calculation of the measure of non-classicality i', it is 
sufficient to consider the above mentioned (/i-range when evaluating the integrals /+ and /_ . 

V. CONCLUSIONS 

We have considered a class of states in an ensemble of two-level atoms, a superposition of two distinct atomic 
coherent states which are called atomic Schrodinger cat states. According to the relative positions of the constituents 
we have defined polar and nonpolar cat states. We have investigated their properties based on the spherical Wigner 
function, which has been proven to be a convenient tool to investigate the quantum interference effects. 

We have shown that nonpolar cat states generally exhibit squeezing, for which we have introduced the measure 
S. The squeezing depends on the separation of the components of the cat and on the number of the atoms the 
cat is consisting of. By solving the master equation of this system embedded in an external environment we have 
determined the characteristic times of decoherence, dissipation and non-classicality of an initial polar cat state. We 
have shown how these depend on the number of the microscopic elements the cat consists of, and on the temperature 
of the environment. Our results show that the decoherence of the polar cat state is surprisingly slow: tdiss/^dec is 
less then a half of an order of magnitude for zero temperature, making these states potentially significant in several 
areas of quantum physics, e.g. experimental studies of decoherence, quantum computing and cryptography. We have 
visualised the process, governed by the interaction with the external environment, using the spherical Wigner function. 
Its transformation in time reflects the characteristics of the behaviour of the atomic subsystem in a suggestive way. 
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